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ABSTRACT 

We explore the extinction properties of the dust in the distant universe 
through the afterglows of high-redshifted GRBs based on the "Drude" model 
which, unlike previous studies, does not require a prior assumption of template 
extinction laws. We select GRB 070802 at z ~ 2.45 (which shows clear evidence 
for the 2175 A extinction bump) and GRB 050904 at z « 6.29, the 2nd most 
distant GRB observed to date. We fit their afterglow spectra to determine the 
extinction of their host galaxies. We find that (1) their extinction curves differ 
substantially from that of the Milky Way, the Small and Large Magellanic Clouds 
(which were widely adopted as template extinction laws in literature); (2) the 
2175 A extinction feature appears to be also present in GRB 050904 at z ~ 6.29; 
and (3) there does not appear to show strong evidence for a dependence of dust 
extinction on redshifts. The inferred extinction curves are closely reproduced in 
terms of a mixture of amorphous silicate and graphite, both of which are ex- 
pected supernova condensates and have been identified in primitive meteorites 
as presolar grains originating from supernovae (which are considered as the main 
source of dust at high-z). 

Subject headings: dust, extinction - galaxies: high-redshift - galaxies: ISM - 
gamma rays: bursts 



1. Introduction 

Dust is present in the high-redshift (z > 2) universe, as evidenced by the reddening 
of background quasars, the depletion of heavy elements in quasar absorption systems, and 
the far infrared (IR) to millimeter (mm) thermal emission of distant quasars. Dust plays a 
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crucial role in the formation and evolution history of stars and galaxies in the early universe. 
The importance of correcting for dust extinction in the universe is now widely recognized. In 
order to reveal the structure and evolution of the early universe, to use Type la supernovae 
(SNe) as standard candles, and to infer the cosmological star formation rate, it is essential 
to correct for the effects of dust extinction. 

Gamma-ray bursts (GRBs), owing to their intense luminosity (emitting up to ~ 10 53 erg), 
allow their detection up to very high redshifts at z > 10 (Lamb & Reichart 2000). Partic- 
ularly, the association of long-duration bursts with massive stars (and therefore with dusty 
regions of high-mass star formation) and the featureless, power law-like spectral shapes of 
their afterglows, make GRBs an excellent probe of the dust at high-redshifts. 

In this Letter we explore the dust extinction of the host galaxies of GRB 070802 at 
z ~ 2.45 and GRB 050904 at z ~ 6.29. We aim at a quantitative examination of the nature 
of the dust in the early universe and attempt to address one of the hotly-debated questions 
in high- 2: astrophysics: do the dust properties evolve as a function of redshift (particularly 
at z > 5 where the dust source may be different)? 



2. Dust Extinction Model 

We characterize the dust extinction properties of GRB hosts with the extinction quantity 
(e.g. Ay r , the rest-frame visual extinction) and the wavelength-dependence of the extinction 
(i.e. A\j 'Ay or A v j Ay if expressed in frequency u, often known as the "extinction curve" 
or "extinction law"). We derive A Vr and A x /A v (or A u /A v ) by fitting the ultraviolet (UV), 
optical, and near-IR afterglow photometry with a dust-reddened power-law model through 



F U = F Q (i//Hz)^exp 



Ay r A(i +Z ) v 



1.086 A 



V r 



(1) 



where F v is the afterglow photometry (with the Galactic foreground extinction corrected), 
P is the intrinsic power-law slope of the afterglow, F is a normalization constant, A( 1+2 )„ is 
the rest-frame extinction, and z is the GRB redshift. 

Unlike previous studies which often assume a template extinction law for A\j Ay (i.e., 
the extinction curves of GRB hosts are assumed to resemble that of the Milky Way [MW], 
the Small Magellanic Cloud [SMC], the Large Magellanic Cloud [LMC], the "Calzetti" atten- 
uation law of starbust galaxies [Calzetti et al. 1994], or the relatively flat "Maiolino" curve of 
AGNs [Maiolino et al. 2001]), we take the "Drude" model proposed in Li et al. (2008a). This 
approach approximates the wavelength-dependence of the extinction by a simple formula 
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consisting of four dimensionless parameters (ci, C2, C3, and C4) 

ci 

(A/0.08) C2 + (0.08/A) C2 + c 3 

233 [1 - cjj (6.88 C2 + 0.145 C2 + c 3 ) - c 4 /4.60] 

(A/0.046) 2 + (0.046/A) 2 + 90 

C A (2) 

(A/0.2175) 2 + (0.2175/A) 2 - 1.95 ' 

where A is in /iin, the first term in the right-hand side represents the far-UV extinction rise, 
the second term and the third term respectively account for the near-IR/ visible extinction 
and the 2175 A extinction bump. 

Compared to models based on template extinction curves, the "Drude" model is pre- 
ferred because (1) it eliminates the need for a prior assumption of template laws - after all, 
there is no reason to assume that the "true" extinction curves of GRB hosts should resemble 
any of those templates, and (2) the analytical formula (eqj2]) on which the "Drude" model is 
based restores the widely-adopted MW, SMC, LMC, "Calzetti", and "Maiolino" templates 
- if the "true" extinction curve of a GRB host happens to resemble a certain template law, 
the "Drude" approach will allow us to restore it (see Li et al. 2008a). 



Ax/Ay = 
+ 
+ 



3. Results 

We apply the "Drude" model to GRB 070802 at z « 2.45 and GRB 050904 at z « 6.29. 
They are selected for the following reasons: (i) they span a wide range of redshifts, from 
the moderately high redshift of z ~ 2.45 (GRB 070802) to the 2nd highest redshift observed 
to date of z » 6.29 (GRB 050904); (ii) the afterglow photometry of GRB 070802 provides 
the most definite evidence for the presence of the 2175 A extinction feature in a GRB host 
galaxy (Kriihler et al. 2008; Eliasdottir et al. 2008); and (iii) the peculiar UKIRT z band 
(Arest ~ 1275 A) flux suppression of the GRB 050904 afterglow at 0.5 days and 1 day after the 
burst (Haislip et al. 2006; Stratta et al. 2007) was interpreted as evidence for an evolution 
of the dust properties at z > 6 (Stratta et al. 2007). 

Using eqs. f|T|2|) and the Levenberg-Marquardt minimization algorithm, we fit the broad- 
band spectral energy distributions (SEDs) of the afterglows of these GRB^] with /3, Ay, c±, 
C2, C3 and C4 allowed to vary as free parametersj§ Therefore, in the SED modeling we have 



1 For GRB 050904 we will consider three different epochs after the burst. 

2 F is not really a free parameter; for a given set of (/?, Ay, c\, C2, C3, C4), F is uniquely determined by 
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six free parameterso It is unfortunate that the number of model parameters (N paxa , = 6) 
exceeds the number of photometry data points iV data for GRB 050904 (iV data = 4 for all three 
epochs; Haislip et al. 2006; Tagliaferri et al. 2005). With iV data = 7, GRB 070802 has a better 
wavelength coverage^ We therefore use xV-^data as a quality measure of the fit. 

In Figure 1 we plot the "Drude" model fit to the afterglow SED of GRB 070802 as well 
as the derived extinction curve. The results for GRB 050904 at three different epochs after 
the burst are shown in Figure 2. We see in these figures that (1) the "Drude" model provides 
excellent fits to the observed SEDs; (2) the derived extinction curves differ substantially from 
the widely-adopted template extinction laws; (3) the 2175 A extinction feature appears to be 
also present in the afterglow spectra of GRB 050904, the 2nd most distant GRB observed to 
date, at epochs of 0.5 days and 1 day after the burst; (4) at an epoch of 3 days after the burst, 
the 2175 A feature appears to be absent in GRB 050904, suggesting that its carrier may have 
been destroyed by the burst jf| and (5) there does not appear to show strong evidence for a 
dependence of dust extinction on redshifts (although the extinction curve does vary from one 
burst to another), as supported by a systematic study of > 20 GRBs at z > 2: the overall 
wavelength dependence of extinction, the steepness of the far-UV extinction rise, and the 
presence and strength of the 2175 A extinction bump, do not appear to show any dependence 
on redshifts (S.L. Liang & A. Li 2008, in preparation). The model parameters are tabulated 
in Table 1. 



the overall flux level. 

3 Admittedly, the models based on template extinction laws have fewer parameters: with the shape of 
the extinction curve fixed, they only need to determine (3 and Ay . The "Drude" approach needs four more 
parameters (i.e. c±, C2, C3 and C4) to describe the wavelength-dependence of the extinction. This is the 
nature of the "Drude" approach; because of this the "Drude" approach is more flexible in revealing the 
"true" extinction curve. 

4 For GRB 070802, we adopt the optical and near-IR photometry of Kriihler et al. (2008) obtained by 
the 7-channel Gamma- Ray Burst Optical and Near-IR Detector (GROND) mounted on the 2.2m ESO/MPI 
Telescope. The ESO VLT spectroscopy of GRB 070802 is in close agreement with the GROND photometry 
(see Fig. 5 of Ehasdottir et al. 2008). 

5 Indeed, one sees in Figure 2 a gradual flattening of the far-UV extinction rise from 0.5 days to 1 day 
and 3 days after burst, as expected from a preferential destruction of small grains responsible for the far-UV 
extinction by the burst (see Perna et al. 2003), that is reflected in Table 2 with a gradual increase (decrease) 
of the cutoff sizes (the power-law size distribution indices). 
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4. Discussion 

In deriving the extinction of GRB hosts, a major problem with the models based on 
template extinction laws is that the wavelength-dependence of the extinction is fixed. For 
a featureless, power-law-like afterglow SED, this often leads to a preference of a SMC- 
type extinction and a small amount of A v (usually < 0.2 mag): obscured by a SMC-type 
extinction (which is roughly a power-law A\ oc A~ L2 ), an intrinsic power-law-like afterglow 
SED remains featureless and becomes a steeper power-law. However, if the dust is "gray" 
(i.e. the extinction A\ only weakly varies with A), the resulting dust-obscured afterglow 
SED will still be a featureless power-law, with the intrinsic power-law exponent unchanged. 
The possible presence of gray extinction has been suggested by a number of authors (e.g. 
see Savaglio et al. 2003, Savaglio & Fall 2004, Stratta et al. 2004, 2005, Chen et al. 2006, Li 
et al. 2008b, Perley et al. 2008). The "Drude" approach allows us to break the degeneracy 
between "gray" extinction and SMC-type extinction. 

Attempts have also been made to fit the afterglow SEDs with the MW, SMC and LMC 
template extinction laws. As shown in Figures 1,2, no acceptable fits are obtained, except 
that the MW model for GRB 070802 and the SMC model for GRB 050904 at an epoch of 

3 days after the burst fit the observed SEDs reasonably well. But even for these two cases, 
the "Drude" approach fits better as can be seen in Figs. 1,2 and indicated by xV^data (see 
Table 1). 

Elasdottir et al. (2008) tried to fit the VLT/FORS2 spectroscopy and the GROND 
photometry with the Fitzpatrick & Massa (1990; hereafter FM) parametrization as well as the 
MW-, LMC-, and SMC-type extinction. They found that satisfactory fits could be achieved 
only if one assumes a cooling break in the intrinsic spectrum, with the FM parametrization 
providing the best fit. However, one should caution that the FM parametrization is only 
valid for A < 2700 A, while the GROND photometry of GRB 070802 extends from ~ 1400 A 
to ~ 6400 A (in the GRB rest-frame). 

The afterglow SEDs of the bursts discussed here all show a flux suppression at A ~ 

4 — 6 iimT 1 and deviate appreciably from a power-law (except for GRB 050904 at 3 days 
after the burst). As shown in Figures 1,2, the flux suppression is closely accounted for in 
terms of dust with a 2175 A bump in its extinction. For GRB 070802, the derived 2175 A 
bump is comparable or even slightly stronger than that of the MW: for MW C4 ~ 0.05 while 
C4 ~ 0.06 for GRB 070802. To validate the suggested detection of the 2175 A extinction 
feature, we have also tried to fit the afterglow SEDs with the "Drude" approach but setting 
c 4 = (i.e. no 2175 A extinction bump). It is found that the fits (with c 4 = 0) are much 
worse, as reflected from the substantially increased xV^data (see Table 1). 



- 6- 



The 2175 A bump, first detected by Steelier (1965), is the strongest spectroscopic inter- 
stellar extinction feature. This feature is seen in extinction curves along lines of sight in the 
MW and LMC0 But it is rarely seen in the afterglow spectra of GRBs. So far, its possi- 
ble detection is only reported in four bursts: GRB 970508 (Stratta et al. 2004), GRB 991216 
(Kann et al. 2006, Vreeswijk et al. 2006), GRB 050802 (Schady et al. 2007), and GRB 070802 
(Fynbo et al. 2007, Kriihler et al. 2008), with the latter showing the clearest presence of the 
2175 A extinction feature in its afterglow spectrum. In addition, Ellison et al. (2006) reported 
the detection of this feature in an intervening absorber at z ~ 1.11 toward GRB 060418. But 
the host galaxy of GRB 060418 at z ~ 1.49 seems to have a SMC-type extinction law. 

The possible detection of the 2175 A extinction feature has been reported for a number 
of low, intermediate, and moderately high redshift systems through (1) the composite ab- 
sorption spectrum of intervening Mgll absorption systems (Malhotra 1997: 0.2 < z < 2.2) or 
radio galaxies (Vernet et al. 2001: z ~ 2.5)0 (2) the individual absorption spectra of interven- 
ing Mgll absorbers (Wang et al. 2004: 1.4 < z < 1.5; Srianand et al. 2008: z ~ 1.3)0 (3) the 
UV SEDs of massive, UV-luminous star-forming galaxies (Noll & Pierini 2005: 2 < z < 2.5; 
Noll et al. 2007: 1 < z < 2.5); and (4) the extinction curves of gravitational lensing galaxies 
(Toft et al. 2000: z « 0.44; Motta et al. 2002: z « 0.83; Wucknitz et al. 2003: z « 0.93; 
Muhoz et al. 2004: z w 0.68). However, Vijh et al. (2003) found that the dust in 906 Lyman 
break galaxies at 2 < z < 4 does not exhibit the 2175 A extinction feature. This is probably 
related to the survival and destruction of the carriers of the 2175 A bump in different physical 
conditions. 

Although the precise nature of the carrier of the 2175 A extinction feature remains un- 
known, it is generally accepted that it arises from small graphitic dust or a cosmic mixture 
of polycyclic aromatic hydrocarbon (PAH) molecules (Li & Draine 2001). In view of the 
detection of presolar graphite dust with a SN origin in primitive meteorites, it is not un- 
reasonable to expect a 2175 A extinction bump for high-z objects since the dust at z > 5 



6 Most SMC extinction curves have no detectable 2175 A bump (Prevot et al. 1984). But there exist 
regional variations in the SMC extinction curve. The SMC sight lines which show no 2175 A bump all pass 
through the SMC Bar regions of active star formation (Prevot et al. 1984; Gordon & Clayton 1998). The 
2175 A bump is seen at least in one line of sight, Sk 143 (AvZ 456), which passes through the SMC wing, a 
region with much weaker star formation (Gordon & Clayton 1998). 

7 But York et al. (2006) found no evidence for the 2175 A bump in the composite absorption spectra of 
809 intervening QSO Mgll absorbers at 1 < z < 1.9. 

8 The 2175 A extinction feature, the 9.7 fim. silicate absorption feature, and the diffuse interstellar bands 
are seen in the damped Lya absorber at z sa 0.524 toward the BL lac object AO 0235+164 (Junkkarinen et 
al. 2004, Kulkarni et al. 2007). 
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is thought to originate from Type II SNe. On the other hand, PAHs have been detected in 
ultraluminous IR galaxies and submm galaxies at z > 2 through their vibrational bands at 
6.2, 7.7, 8.6 and 11.3 /im (see Lutz et al. 2005, Yan et al. 2005). PAHs were also seen in the 
Cloverleaf lensed QSO at z ~ 2.56 (Lutz et al. 2007). If PAHs are indeed responsible for the 
2175 A extinction, it would not be surprising to see this feature in high-z galaxies. 

Finally, we fit the inferred extinction curves using a mixture of spherical amorphous 
silicate and graphite dust each with an exponential-cutoff power-law size distribution (e.g. 
see Kim et al. 1994) 

/'"■max 

A x /A v = A sil Cf t (a,A)a- Qsil exp (-a/a CjSil ) da 

/^max 

+A gra / C c g ™(a, A) a~ a ^ exp (-a/a c , gra ) da, (3) 

where the lower (upper) cutoff size a m \ n (a ma x) is taken to be 50 A (1 /im) for both silicate and 
graphite dust; the power-law indices a s n, a gra and the exponential-cutoff sizes a CjS ii and a C)gra 
are treated as free parameters; A s \\ and A gT& are related to the abundance of each species; 
and Cgxt (CSt) is the extinction cross section of silicate (graphite) dust. As shown in Figure 
lc and Figure 3, the silicate-graphite model closely reproduces the inferred extinction curves 
for both GRBs, including the 2175 A extinction bump (see Table 2 for the size parameters). 
The major mismatch occurs at A ~ 7/mi~ 1 which is probably due to the sudden rise of 
the silicate electronic absorption (see Kim & Martin 1995). We note that both silicate and 
graphite are expected SN condensates (Todini & Ferrara 2001, Nozawa et al. 2003). They 
have been identified as presolar grains in primitive meteorites originating from supernovae 
which are considered as the main source of dust at z > 5 (see Dwek et al. 2007). 

By fitting the afterglow SEDs of GRB 050904 (z « 6.29) with the extinction curve 
inferred for the distant BAL QSO at z ~ 6.2 (which displays a plateau at A -1 ~ 3.3 — 
5.9 /im -1 , Maiolino et al. 2004), Stratta et al. (2007) argued that the dust properties may 
evolve beyond z > 5. This seems to be supported by that the dust at z > 5 is probably 
produced by Type II SNe while in the local universe AGB stars are a major source of dust. 
However, this study together with a preliminary analysis of > 20 GRBs at z > 2 based on 
the "Drude" approach does not indicate any dependence of the dust extinction on redshift. 
A more thorough and systematic study of the dust extinction and IR emission properties of 
high- z GRBs is in progress and will be used to further explore whether the dust properties 
vary as a function of redshift. 

We thank V.P. Kulkarni, S. Savaglio, G. Stratta, D. Watson and D.M. Wei for very 
helpful comments and suggestions. We are supported in part by a NASA/Swift Theory 
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Table 1 : Parameters for fitting the afterglow SEDs with the "Drude" model and the MW, LMC and SMC 
template extinction laws. 
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Table 2. Dust size distributions for the extinction curves derived from the "Drude" model 

and a mixture of silicate and graphite grains 
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Fig. 1— Left panel (a): Fitting the SED of the afterglow of GRB 070802 (z « 2.45) with the 
"Drude" approach (red) and the MW (magenta), LMC (blue) and SMC (green) templates 
for the GRB host extinction curve. Middle panel (b): Comparison of the MW (magenta), 
LMC (blue), and SMC (green) extinction laws with that derived from the Drude approach 
(red). Right panel (c): Fitting the derived extinction curve (red solid line and black filled 
circles) with a mixture of amorphous silicate (cyan dotted line) and graphite dust (green 
dashed line). The blue solid line plots the resulting model extinction curve. 
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Fig. 2. — Same as Figure la,b but for GRB 050904 (Haislip et al. 2006; Tagliaferri et al. 
2005) at three different epochs after burst. 
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Fig. 3. — same as Figure lc but for GRB 050904 at three different epochs after burst. 



